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ABSTRACT 

This letter presents a theoretical derivation of an estimate for a radio source 
jet kinetic luminosity. The expression yields jet powers that are quantitatively 
similar to a more sophisticated empirical relation published by the Willott, Blun¬ 
dell and Rawlings at Oxford. The formula allows one to estimate the jet kinetic 
luminosity from the measurement of the optically thin radio lobe emission in 
quasars and radio galaxies. Motivated by recent X-ray observation, the deriva¬ 
tion assumes that most of the energy in the lobes is in plasma thermal energy 
with a negligible contribution from magnetic energy (not equipartition). The 
close agreement of the two independent expressions makes the veracity of these 
estimates seem very plausible. 

Subject headings: quasars, jets , radio galaxiesgalaxies:active — galaxies:jets - 
quasars: general 

1. Introduction 

The purpose of this letter is to discuss estimates of the power transported by the radio 
jets in quasars and radio galaxies. An accurate estimate of the jet power is of fundamental 
physical interest, since it can be used to quantify the power emerging from the central engine 
of the radio source. In actuality, the radio luminosity is merely an indirect measure of the 
energy transported through the jets from the central engine that is not readily interpretable. 
Most of the energy flux is in mechanical form (kinetic luminosity) - the particles and fields 
necessary to produce the synchrotron luminosity that is detected in the radio lobes. The 
radiation losses, manifested as radio emission from the jet, are merely the waste energy of 
this kinematic flow. 


Surprisingly, the most difficult methods of estimating jet power rely on observations of 
the jets themselves. Due to significant Doppler enhancement in relativistic jets, the syn¬ 
chrotron radio emission is a poor indicator of intrinsic jet power. For example, Cygnus A 
has extremely powerful radio lobes and faint radio jets. Most of the energy in the jets is not 
radiated away, but is transported to the lobes in the classical FRII double lobe morphology. 
Even the inclusion of observations of high energy emission such as optical or X-ray (inverse 
Compton) in one’s analysis of jet energetics does not tightly constrain the bulk jet flow. If 
the resolution is poor at high frequency (as is often the case) then one can not necessarily 
associate the plasma emitting the high frequency photons with the radio emitting plasma. 
If one has high resolution images then the high frequency emission can be detected in en¬ 
hancement regions or knots in the jets. One can use this information to get an estimate of 
the plasma conditions within the dissipative knot, but this do not necessarily constrain the 
plasma state in the bulk of the jet. Furthermore, there are still ambiguities with the Doppler 
factor that affect the estimates quite dramatically. 

The Doppler enhancement of relativistic flows in jets is a crucial parameter since the 
total luminosity of an unresolved jet scales as the Doppler factor to the fourth power and to 
the third power for a resolved cylindrical jet (Lind and Blandford 1985). This is the reason 
why the implementation of 5 GHz flux densities, as is common in studies of radio loudness 
of large quasar samples, is a poor indicator of the true intrinsic kinetic luminosity of the 
jets. More specifically, the majority of core dominated blazar-like quasars have incredibly 
strong 5 GHz flux densities from emission on the subkiloparsec scale, yet they have weak or 
moderate radio lobe emission Punsly (1995). This is interpreted as the jet being of modest 
kinetic luminosity (at most) because there is not a large amount of hot plasma and gas that 
has been transported through the jets to the radio lobes. The 5 GHz flux only represents the 
dissipation in the jet itself and it has been extremely Doppler boosted. An estimate of kinetic 
luminosity based on 5 GHz flux density can be off by four or more orders of magnitude for 
a core dominated blazar. 

A far better way to estimate the kinetic luminosity from a jet is to study the isotropic 
properties of the material ejected from the ends of the jets in the radio lobes. The radi¬ 
ation from the lobe material is generally considered to be of low enough bulk velocity so 
that Doppler enhancement is not much of an issue. The basic idea is that lobe expansion 
is dictated by the internal dynamics of the lobes and the physical state of the enveloping 
extragalactic gas. X-ray observations can indicate a bremsstrahlung spectrum of the sur¬ 
rounding gas that can be used to find the pressure of the extragalactic medium. X-rays also 
provide information on the working surfaces at the end of the lobes, ’’the hot spots.” One 
can associate the X-ray emission as inverse Compton radiation from the hot spots and the 
radio luminosity is the synchrotron emission from the hot spots. This constrains the plasma 
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state within the luminous hot spots. However, most of the energy stored in the lobes is in 
the large diffuse regions of radio emission that constitutes the majority of the large volume 
of the radio lobes. It is the enormous volume of synchrotron emitting plasma within the 
lobes (~ 10 4 — 10 5 kpc 3 ) that is the most direct indicator that the jets must be supplying 
huge quantities of hot plasma and magnetic field energy to the lobes. One can also use the 
curvature of the lobe synchrotron radio spectra to estimate parameters in the diffuse lobe gas 
- this is known as spectral ageing. Of course, all of these plasma state estimations are most 
accurate when applied to situations in which one has deep X-ray and radio data of a relaxed 
classical double lobe structure. This only occurs in a few instances, so such detailed analysis 
are not compatible with large sample studies. Motivated by these limitations, this letter 
presents two techniques for estimating jet energy based on partial information on the lobe 
parameters. The two methods involve different assumptions and have different ambiguities. 

The most sophisticated calculation of the jet kinetic luminosity incorporates deviations 
from the minimum energy estimates in a multiplicative factor / that represents the small 
departures from minimum energy, geometric effects, filling factors, protonic contributions 
and low frequency cutoff (see Willott et al (1999) for details). The quantity, /, is argued 
to constrained to be between 1 and 20. In Blundell and Rawlings (2000), it was further 
determined that / is most likely in the range of 10 to 20. Therefore we choose / = 15 in 
order to convert 151 MHz flux densities, F\^\ , to estimates of kinetic luminosity, Q 151 , using 
equation (12) and figure 7 of Willott et al (1999), 

Q 151 ~ 1.1 X 10 45 [(1 + z) 1+a Z 2 F 151 ] 1 ergs/sec , (1-1) 

Z « 3.31 - 3.65 x 

([(1 + zf - 0.203(1 + zf + 0.749(1 + zf + 0.444(1 + z) + 0.205] "°' 125 ) . (1-2) 

The quantity F 151 is the optically thin flux density from the lobes (i.e., no contribution from 
Doppler boosted jets or radio cores) measured at 151 MHz in Jy. The flux density spectral 
index is defined as F u ~ v~ a . We have assumed a cosmology with F[ o=70 km/s/Mpc, 
Da = 0.7 and Q m = 0.3. The expression for Z is from Pen (1999). 

In the following, a new estimate of the jet kinetic luminosity is derived that is motivated 
by the wealth of X-ray data on radio lobes that has been published since Willott et al 
(1999); Blundell and Rawlings (2000). Both the current manuscript and the Willott et al 
(1999) derivations rest on the basic premise that Q = U/T + radiation losses, where U 
is the energy stored in the lobes and T is the elapsed time. In Willott et al (1999), U is 
found by assuming the lobes are near equipartition and there is uncertainty in the energy 
from protonic components and low frequency cutoffs which are incorporated in the empirical 
factor / discussed above. Conversely, motivated by the new X-ray data described below, U is 
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computed theoretically in the limit that the lobes are very far from equipartition. In Willott 
et al (1999), T is determined by an empirical estimator for the age of the radio source based 
on their lengths and head advance speeds. Conversely, in this treatment T is computed from 
spectral ageing. In spite of the fact that U and T are determined from completely different 
methods and assumptions, the expressions for Q that are found in (1.1) and (3.9) yield 
similar values (to within a factor of 2) for the jet luminosity. This close agreement lends 
credence to the claim that these formulae are robust estimators of jet kinetic luminosity. 


2. Motivation: X-ray Observations 

The minimum energy condition in the lobes seems to be in conflict with the X-ray 
data on the surrounding extragalactic gas. This was first noted for Cygnus A (see Punsly 
(2001) and references therein) and later for 3C 388 in Leahy and Gizani (2001) as well as 
for a large sample of FRII radio sources in Hardcastle and Worral (2000) based on ROSAT 
data. These studies concluded that typically the pressure in the external gas greatly exceeds 
(by at least an order of magnitude) the lobe pressure associated with the minimum energy 
assumption. The general picture that seems to be emerging from the X-ray data of ROSAT, 
ASCA, XMM and Chandra is that the hot spot energies seem to slightly exceed the minimum 
energy requirement based on inverse Compton calculations of X-rays from the hot spots (see 
for example Wilson et al (2001); Brunetti et al (2002b)), but the lobes themselves are far 
from equipartition. In Leahy and Gizani (2001), it was deduced that the departure from 
equipartion was most likely the consequence of a low energy population of positrons and 
electrons that is not an extension of the power law distribution responsible for the radio 
emission, but a low energy excess. The other possibility is protonic matter which would also 
drastically increase the energy content of the lobes over the equipartition estimates. In either 
case, the held energy is only a few percent of the particle energy in the lobes. For example, 
high resolution Chandra data was used to model the X-rays as inverse Compton emission in 
the lobes in the radio galaxy 3C 219 (Brunetti, G. et al 2002a). It was concluded that the 
particle energy exceeded the magnetic energy by a factor of 60 in the radio lobes. Similarly, 
the FRII radio galaxy 3C 452 was studied with Chandra in Isobe et al (2002) who estimate 
that the energy density in the particles is 27 times that of the magnetic held in the lobes. 
For most FR II sources, typical magnetic held strengths in the lobes that are estimated from 
X-ray data and pressure balance are only a third to a fifth of the minimum energy value 
(Hardcastle and Worral 2000; Leahy and Gizani 2001). 

These X-ray observations yield valuable information on the energy content of the lobes, 
but are disjoint from the spectral ageing estimates of the lobe advance speeds. In the 
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following estimation of jet kinetic energy, it is assumed based on the X-ray data presented 
above that the energy content of the lobes is purely in particle form to first order (accurate 
to a few percent). Yet, the notion that spectral ageing provides an estimate of lobe age is 
retained. By choosing a subequipartion held strength in the lobes, spectral ageing estimates 
are found to be longer than the corresponding minimum energy estimates (Alexander and 
Pooley 1996). Similarly, by setting t sep = t syn , the subequipartion fields yield lower lobe 
separation velocities if this spectral ageing argument is viable. Thus, the problem of the large 
lobe advanced speeds in the minimum energy assumption is remedied by this modification 
(Alexander and Pooley 1996). This method of computing jet kinetic luminosity is the lowest 
order improvement to the minimum energy estimate and was implemented in Punsly (2001) 
to study Cygnus A. 


3. Particle Energy Dominated Lobes 

Motivated by the X-ray observations, we proceed to compute the jet power based on 
the limit that all of the lobe energy is in the hot particles. We also assume that the time to 
convert a jet energy flux to the stored lobe energy is the time that it has taken the lobes to 
propagate from the central engine to their current separation, t sep . 


3.1. Spectral Ageing 


Spectral ageing within the radio lobes is often used to determine the lobe plasma age. 
The results are predicated on the assumption that the lobe plasma is primarily back flowing 
plasma in the sense that jet plasma is deflected backward at the working surfaces in the hot 
spots to form the lobe plasma. By studying the curvature of the radio spectra at different 
points within the lobes, one can in principle (if there is no reheating or re-injection of the 
plasma) determine the gradient in the high energy cutoff of the electron distribution due 
to synchrotron cooling and hence the plasma age. The age of the lobe plasma closest to 
the central engine should be the oldest plasma. Thus, one has an estimate of lobe age 
and therefore the lobe advance speed. Defining the spectral break frequency as u b , the 
synchrotron lifetime is expressed in cgs units as (Liu et al 1992) 


t 


syn 


1.58 x 10 12 5 -§z/~5 . 

b 




(3-1) 
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3.2. The Energy Contained Within the Synchrotron Emitting Plasma 


Consider a power law distribution of energetic particles (probably electrons and positrons) 
expressed in terms of the thermal Lorentz factor, 7 , for a uniform source in a volume, V. 
The total number of particles contributing to the synchrotron radiation in the frequency 
interval u 1 < v < v 2 is: 

N r = N 0 V /Vd 7 . (3-2) 

Jji 

The minimum and maximum Lorentz factors in the expression above are related to lower and 
upper cutoff frequencies in the synchrotron spectrum and u 21 respectively by (Ginzburg 
1979) 
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where u B = (eB )/(27rm e c) is the cyclotron frequency and note that 


(3-3) 


yM) = 2-2, y 2 ( n ) = 0 . 10 , if n = 2.5 : and y,(n) = 2.7, y a (n) 


0.18, if n = 3.0 .(3-4) 


The synchrotron spectral luminosity of the plasma, L(v), is a function of both the particle 
distribution in momentum space and the magnetic held strength. Integrating the synchrotron 
power formula over the particle distribution yields (Ginzburg 1979) 


where 


U e 


2 x 10 U B~I 
a(n)(n — 2 ) 


/ \ _1 r . \ n— 1 

L ("i 2 


1 

1 

( hi 

n—1 “ 

[ n ) v i \ 2 

L 


(«) z/ 2 / 


(3-5) 


a(n ) = 


( 2 “f 1 v / 3 ) r (2!gi) r (2s±i») r (qi) 
8 v/?(n + l)r (77 


(3-6) 


3.3. Estimating the Jet Power 

Set t sep equal to the synchrotron ageing timescale t syn , associated with the spectral break 
in the flux density, F u , of the lobe plasma closest to the quasar (the emission just above the 
spectral break is from the lowest energy electrons that have synchrotron radiated away their 
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energy and hence the oldest subpopulation of charges that have experienced synchrotron 
decay in the lobes). By combining t syn from (3.1) with the expression for the plasma energy, 
(3.5), one obtains an estimate for the energy stored in the lobes as a function of spectral 
luminosity, L(u), in the limit of particle energy dominance, 


U P 


rMM 1 / 2 W 1/2 

7.9 (n — 2 )a(n) 


[yi(n). 


^ — 1 
2 t 


syn 


(3-7) 


Evaluation of the formula above requires numerous characteristic frequencies that need 
to be determined. Since expressions that are applicable to sparse data are desired for eval¬ 
uating large samples, we choose a common set of “typical” parameters for an FRII radio 
source. First of all, determining requires high resolution maps of the lobes at a variety 
of frequencies. This data has been obtained for only a limited number of bright sources. 
The largest sample of these detailed observations is from Liu et al (1992). The average rest 
frame break frequency from the sample of Liu et al (1992) is zy, = 8.9 ± 7.0 GHz. Secondly, 
in order to estimate the minimum synchrotron frequency we note that from Braude et al 
(1969) (even though the measurements arc likely to be extremely inaccurate) it is clear that 
many FRII sources are very strong emitters down to frequencies at least as low as 12.6 MHz. 
Thus, we pick = 10 MHz in the quasar rest frame. Finally, in order to approximate the 
total radio luminosity, L = f L{y) dv (including the significant contribution at frequencies 
above the spectral break), with a single spectral index, a value of z/ 2 = 100 GHz is chosen. 
Inserting these “typical” frequency values into (3.7), one obtains a simple estimator of lobe 
power in the limit of particle dominance, and noting that at the spectral break frequency, 
t t 

L syn ~ b sep 


Q 


u. 


+ L 


v sep 


n — 1 

[ yi (n)]—( 15.1) a 
(n — 2 )a(n) 


10 42 (1 + z) 1+a Z 2 Fi 5 i ergs/sec + L , 


(3-8) 


where the spectral index a — (n — l)/2 has been introduced ( L[y ) ~ u~ a ). It should be 
noted that the estimates above are very conservative. The existence of a substantial proton 
component to the lobe gas or an extension of the low frequency portion of the electron 
spectrum would increase the energy flux estimates significantly. Observations suggest that 
a ~ 1 is a good fiducial value for the expression (3.8)(Kellcrmann et al 1969), 


Q par « 5.7 x 10 44 (1 + z) 1+a Z 2 Fi 5 i ergs/sec , a^l. (3-9) 


4. Conclusion 

An independent formula for the jet kinetic luminosity estimator in (1.1) is derived 
in (3.9) that was motivated by different physical assumptions. The two estimates agree to 



within a factor of 2. This lends credence to the idea that (1.1) and (3.9) are robust estimators 
of jet kinetic luminosity when the optically thin extended emission is measured in a deep 
radio map. The main result of this paper is that (1.1), although very ambitious in its intent, 
is likely to be correct to within a factor of a few even if some of the assumptions in its 
derivation are inaccurate. 
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